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Abstract 

A model predictive control (MPC) scheme for a permanent- 
magnet synchronous motor (PMSM) is presented. The torque 
controller optimizes a quadratic cost consisting of control error 
and machine losses repeatedly, accounting the voltage and 
current limitations. The scheme extensively relies on optimiza- 
tion, to meet the runtime limitation, a suboptimal algorithm 
based on differential flatness, continuous parameterization and 
Unear programming is introduced. 

The multivariable controller exploits cross-coupling effects 
in the long-range constrained predictive control strategy. The 
optimization results in fast and smooth torque dynamics while 
inherently using field-weakening to improve the power ef- 
ficiency and the current dynamics in high speed operation. 
As distinctive MPC feature, constraint handling is improved, 
instead of just saturating the control input, field weakening 
is applied dynamically to bypass the voltage limitation. The 
performance of the scheme is demonstrated by experimental 
and numerical results. 

I. Introduction 

The efforts of implementing predictive controllers in elec- 
trical drives aim at replacing the classical cascaded field- 
oriented control structure with PI controllers. The machine can 
be better exploited by improved control behavior, the system 
variables are optimized. In this contribution, the conventional 
torque and current control structure of two separate controllers 
is changed to multi-input multi-output (MIMO) control. By 
transformation into the field-oriented frame, torque generation 
is decoupled from flux variation, however, the current dynam- 
ics are still strongly coupled, therefore a MIMO controller 
is advantageous. The aimed improvements are better current 
and voltage constraint handling by exploiting cross-coupling 
between the orthogonal components, and better power effi- 
ciency and dynamics by optimally adjusting the currents in 
both dynamic and steady-state operation. 

The major obstacle in implementing predictive control 
schemes is the limited computational power, inherited by 
the high sampling rates. The most widespread schemes trade 
computational feasibility against compromises in the prob- 
lem formulation, for instance, generalized predictive control 



(GPC) has a high prediction horizon but is unconstrained, 
whereas predictive torque control (PTC) is constrained but 
so far only reaches 2 steps of prediction fill. To obtain the 
advantages claimed by classical MPC on constrained MIMO 
systems, both, inclusion of constraints and a high prediction 
horizon are required. Only few schemes so far satisfy these 
two requirements 0. Using continuous control in the field- 
oriented frame, the analytical problem description enables us- 
ing efficient optimization algorithms to maximize the obtained 
information for a given computational power 

The online solution of the linearly constrained linear- 
quadratic problem, typical for MPC, requires quadratic pro- 
gramming (QP) algorithms, which are, however, computation- 
ally too expensive for drive systems. A recent development 
is the use of explicit MPC, where an offline solution is 
computed and stored as look-up table in the real-time con- 
troller ||3], The scheme reaches 5 prediction steps with 
constraints. Recently an online algorithm based on a fast- 
gradient algorithm was proposed for reference tracking control 
of a grid-connected inverter ||5]. While the runtime results 
are convincing, experiments are pending. The advantage of 
online optimization in MPC is the possibility to manipulate or 
adapt parameters, which results in simplified commissioning, 
an aspect that is becoming more and more important in the 
recent developments [6j. 

This implementation entitled 'classical MPC relies exten- 
sively on online-optimization. It is embedded in a flatness- 
based predictive control scheme Q, JH], where trajectory 
optimization and prediction are decoupled, and the trajectory 
generation can be simplified ID. A suboptimal trajectory 
generation algorithm [lO] is applied which is based on a 
continuous approach, the variables are not discretized but 
represented as a polynomial with undetermined coefficients. 
The cost function is linearized in the unconstrained optimum 
in order to use a linear programming (LP) solver, which is 
amongst the simplest and fastest numerical optimizers. As re- 
sult, a (suboptimal) prediction of 2 ms with current and voltage 
constraints is obtained at 8 kHz sampling rate. This way, by 
solving the computational problem without compromising the 
MPC formulation, the merits of this controller can be studied. 

The initial results have been presented in ifTTl . 

11. Problem Statement 

A. Machine Model 

In a fist step the machine model is linearized. This will 
result in a linear-quadratic optimization problem with linear 
constraints, a very common problem in optimal control sys- 
tems 1(121 that is also simpler to solve in real-time. Assuming 



that the rotor speed does not change too much over the 
optimization horizon T, 



wmW= const. Vie[0,r], (1) 



the PMSM model and the voltage equations become linear 
The assumption is justified if the current control loop is faster 
than the speed control loop. The electrical subsystem of the 
machine, consisting of the quadrature and direct currents iq 
and id (peak values), is given as 

-Rid + ripUJMLqiq + Ud, (2) 



Ld-id 

L ^^ 
^^dt*« 



TM -UpKiq. 



-Riq ~ UpUiMLdid - ripUiMK + Uq, (3) 

(4) 



The nomenclature is shown in table II in the appendix. It is 
noted that the reluctance torque t^j ~ ^np{Ld — Lq)idiq is 
neglected, as this term is very small compared to the electro- 
magnetic torque in surface-mounted PMSMs or in machines 
with small saliency. Furthermore, it would render the model 
nonlinear, requiring nonUnear optimization methods llT3l . 

B. Optimization Goals and Cost Functional 

The formulation of a suitable cost functional is a key point 
in predictive control, as it is the only tuning possibility of 
the control scheme. The optimization is aiming at minimizing 
the control error for good dynamical performance as well as 
machine losses for better efficiency. Both goals are included 
in the cost functional. By choosing the cost functional and 
weights well, it is possible to find a good trade-off between 
both goals during transients, and eventually to fulfill both goals 
in steady-state. The cost functional for the predictive torque 
controller is 

J = I {Pctrl{t)+WL- Plossit)) dt + T-P,trl{T) (5) 

Jo 

which trades off the squared control error from the constant 
torque reference t'^j 



Pctrlit) = (tm - tIi)^ 



(6) 



with machine losses 



Plossit) - -Ritl + tl) + -npWMfcFeC*' + 

3 3 

= 2^(*d + ''q) + -jJ^P^MkpeiiLdid + K)"^ + iLqiq)'^). 

(7) 

The first term in Pioss represents copper losses, and the 
second term represents the iron losses consisting of hysteresis 
losses. The last term in (|5]l is the end-weight of the control 
error, it serves to improve convergence within finite time, 
especially when constraints are active. Its effect is shown in 
the results section. Eddy current losses are negligible on the 
tested machines, however, they could be included using the 
model presented in llT4l . The machine losses can be reduced 
by field-weakening, where a trade-off between copper and iron 



losses is found ifTsl . Imposing a negative direct current id, the 
flux magnitude in the stator is reduced while the copper losses 
increase. As the iron loss constant kpe is not part of standard 
motor parameters, it has to be determined experimentally llT6l . 
The iron losses are quite considerable, for instance, the applied 
2.64 kW machine has about 43 W copper losses and 200 W 
iron losses in rated operation, therefore the parameter kpe can 
be determined using only standard equipment. Inverter losses 
remain undiscussed as the switching frequency is arbitrary 
fixed. 

There are two tuning parameters, Wl and T. The weight 
Wl was set 0.05, the value was determined heuristically. The 
optimization horizon is set T = 2 ms such that the cost func- 
tional includes the complete setpoint change. It is important 
that the optimization horizon is high enough, otherwise the 
open-loop and closed-loop trajectories differ and the behavior 
is strongly suboptimal. This is illustrated in Fig. [T] (left). 
The open-loop trajectory, a continuous trajectory optimal with 
respect to J, is generated and a first part is applied to the 
system. At the next sampling step (indicated by a circle), a 
new open-loop trajectory is generated and again its first part 
is applied to the system. The resulting closed-loop trajectory, 
which is the sequence of the discrete points at each sampling 
interval, differs significantly from the open-loop trajectories, 
especially if the horizon is too small. Then, the closed-loop 
trajectories simply do not fit the cost functional anymore and 
are suboptimal. For a horizon higher than required for the 
setpoint change, as shown in Fig. [T] (right), the difference 
between open- and closed-loop trajectories becomes smaller. 
Ideally the points should be superposed with the continuous 
open-loop trajectories, then, the closed-loop trajectories can be 
assumed optimal regarding the open-loop cost functional J. 

C. Current and Voltage Constraints 

The most important nonlinearities of a PMSM, in terms 
of control, are the voltage and current limitations. The current 
constraints prevent overheating of the machine, and the voltage 
is limited by the maximum output voltage of the voltage 
source inverter The voltage constraints limit rotor speed 
as well as current dynamics in high-speed operation. Both 
constraints are approximated by affine equations, in order to 
be computationally efficiently treated. 

The current range for the direct current id is limited to 
jmm < < 0. Only negative values of id are desirable, 
as they improve power efficiency and reduce the induced 
voltage by weakening the flux magnitude in the stator ID, 
flSl . The lowest value i™™ is the optimum value at rated 



speed (g-Pi 



0) and is given as 

LdK 



Lj + RinpUJMNkpe 



(8) 



which is independent of quadrature current iq if the reluctance 
torque and the saliency are neglected, what is acceptable for 
surface-mounted PMSMs ifT?! . The value is doubled to enable 
further field-weakening to improve dynamics in high speed, 
an effect described in the results section. For the quadrature 




Fig. 1. Exemplary torque setpoint change to describe open- and closed-loop trajectories in long-range predictive control. Left: small horizon, right: high 
horizon. Circles: starting points for each trajectory planning iteration. 



current iq, the largest possible range of values should be avail- 
able. The resulting linear constraints, shown in Fig. |2] almost 
completely fill the current region of interest. A representation 
with affine inequality constraints is thus acceptable. 

The approximation of the voltage constraints is a bit more 
difficult. The g-axis should not be restricted, as the induced 
voltage is aligned to it and is the largest value that will appear. 
A steady-state analysis of the system equations (2), (3) shows 
that a rectangular voltage area results 

D -Tiiin _ r rnax -max ^ ^ ^ T max -max 
Rl^ -UpLqiO^j <Ud< UpLqUJ^j 1^ , 

(9) 
(10) 

This rectangle (dark grey on Fig. |2) is expanded such that 
the outer circle of the voltage limitation is hit (light grey 
on Fig. lU. During dynamical transients, the voltage vector 
points to one of the outer corners, subsequently touching the 
outer limiting circle. Therefore, a linear approximation of the 
voltage limits as a rectangle by the presented method, as shown 
on Fig. |2l does not limit the steady-state operational range and 
only marginally affects dynamics. A less restrictive method is 
presented in |i5j, where a time-varying constraint in form of a 
hexagon in stator frame is proposed. The limiting circle can 
be expanded to a hexagon by using overmodulation techniques 
ifTTl . While this is possible too with the underlying predictive 
control algorithm, the method in the (d, (7)-frame is chosen for 
simplicity and to prevent possible current ripples. 

III. Optimal Control Algorithm 

To study real-time applicability of the presented scheme, 
first, the highest possible amount of optimization parameters 
is determined. The fastest optimizer with constraints is the 
widely known linear programming (LP) method. Table I shows 
some worst-case computational results of LP (simplex method 
from ifTSl ) as function of the number of free parameters (CPU: 
1.4 GHz industrial PC). More parameters lead to a higher 
number of iterations which are also more complex; the worst- 
case number of iterations is the number of parameters plus the 
number of constraints (T2|. As in the underlying application. 




Fig. 2. Affine approximation of the current and voltage constraints. 
Circle: feasible set of current and voltage vectors, grey: feasible set after 
approximation of the constraints. 

the constraints are decoupled, however, this worst-case is 
not to be expected. The maximum runtime is given by the 
sampling rate minus latency of input/output, therefore at 8 
kHz sampling rate, it must be less than about 110 /is. Thus, at 
best, 12 parameters can be optimized if a LP method is used. 
Runtime of the predictive controller is further discussed at the 
end of the section. 

TABLE I 

Runtime of a linear program for some worst-case problems on 
A 1.4 GHz CPU 



Parameters 


Constraints 


Iterations 


Runtime [/^s] 


20 


44 


67 


769 


12 


28 


34 


165 


8 


20 


10 


35 



A. Trajectory Generation 

The trajectory generation algorithm presented in [ITOl . a 
development related to flatness-based methods |j9|, is chosen. It 
can optimize a quadratical cost function with linear constraints. 
As major differences to standard algorithms, it is applying 
a continuous parameterization instead discretization, and the 
computationally efficient linear programming solver is used 
instead of quadratic programming or iterative gradient search. 



Even though LP is used, it is still quadratic optimization; the 
unconstrained solution to the quadratic cost is calculated first, 
then, constraints are included with the LP solver 

The trajectories for the current are defined as degree n 
power series with undetermined coefficients a^ . 



k=0 



k=0 



rpk ' 



(11) 



This definition reduces the dimensionality of the generated 
trajectories rather than their length and is referred to as Ritz 
parameterization ||121. It is an alternative parameterization to 
the typical Euler discretization. The first coefficients ado and 
ttgo are the initial conditions, and the remaining 6 coefficients 
are determined by optimization. A high prediction horizon is 
obtained for a relatively small number of parameters. Due 
to the analyzed computational limitations, n = 3 is chosen 
as polynomial degree. Fig. [3] illustrates the computational 
advantage, with 3 parameters at 8 kHz sampling rate, using 
a discrete description, the prediction horizon is 0.375 ms, 
but with a degree 3 polynomial, a well-conditioned setpoint 
change can be described over the desired prediction horizon 
of 2 ms. A more complex trajectory may not be expected for 
the underlying application. 



L25 




0,75 1 1,25 
time [ms] 

Fig. 3. A trajectory described by 3 free parameters. Black: polynomial- 
based curve, grey: discrete-time curve. Discrete-time horizon is 0.375 ms, 
continuous polynomial-based trajectory is well-conditioned at the desired 
horizon length 2 ms. 



The corresponding voltages Udq{t) are computed by alge- 
braic differentiation of (fTTT i and by solving for the model 
equations (|2|l and ([3]l, this is also called the flatness-based 
approach |fT9l . This way, the voltages do not need to be 
represented by additional parameters, and equality constraints 
are avoided in the optimization. The variables idq, j[idq and 
Udq are substituted in the cost functional J ^ hy the found 
functionals. 

The cost functional J is then a quadratic function of 
the unknown parameters a and of the parameters which 
were assumed constant, namely the motor parameters, the 
measured currents, the speed ujm and the torque reference 
T^j. Defining the vector of undetermined coefficients as a 



{adi,ad2,ad3,aqi,o:q2,aq3)'^, it is written as 
J = OL^Qa + q^a + go- 



(12) 



Because of the parameterization with a polynomial basis, 
convexity must be discussed [[9], in (|5]l the weight matrix was 
1^2x2 vvhereas now it is extended to Q G K^^^. Convexity 
guarantees the existence of a unique global minimum, and 
convergence of the optimization algorithm in finite time. The 
proof that convexity is maintained with this transformation 
is given in ifTOl . Graphically, for a simple example with two 
parameters, J can be represented as in Fig. |4](left). 

The inequality constraints are also parameterized with the 
polynomial. Exact parameterization requires linear matrix in- 
equality (LMI) methods, which are however too involved for 
this application. A simpler way is to sample the trajectories 
for idq{t) and Udq{t) at an interval ^. In fTOl it is proven that 
if an additional interlay is added, the transformation guaran- 
tees maintenance of the original constraints. For instance, a 
constraint id{t) < OVt G [0,T] is parameterized as 



«d(0) < 0, 



id{k-) - AidiO) <0fc = L.3, 
n 



(13) 
(14) 



where A = 0.064 is the respective interlay constant for n = 3 
(see ifTol ). The first constraint on the initial condition does 
not need to be included, and the remaining 3 conditions are 
affine functions of a, but not of t, such that they can directly 
be included in linear-quadratic optimization in the parameter 
space a. 

As J is convex, the unconstrained global optimum ckq is 
found algebraically by solving first-order necessary conditions. 



1 



Q 



Then, by an affine coordinate transformation 
(3 = A{a. — ckq), 



(15) 



(16) 



the problem can be reformulated as least-distance problem, i.e. 
a quadratical cost describing the distance to the unconstrained 
optimum. The linear transformation includes a shift of the 
origin of the coordinates, a coordinate rotation as well as a 
coordinate scaling, and is found with 



A^A = Q, 



(17) 



which can be solved with the Cholesky decomposition ~ 
cholesky((5-'^). As result, the cost functional looks much 
simpler and reduces to a sum of squares 



J = f3^(3, 



(18) 



and can be represented as in Fig. |4] (middle). The uncon- 
strained optimum is thereby /3g = 0. This constrained least- 
distance problem is already simpler to solve than the original 
problem. In the next step, the least-distance problem is lin- 
earized around the unconstrained optimum, see Fig. E] (right). 



The squares in the cost function are replaced by absolute 
values 

J = /3^/3 = ^/32«^|/?,| = J'. (19) 

i i 

In the LP standard form, furthermore, only positive param- 
eters are possible, therefore the variables are replaced by 
Pi = Pip - Pin, with Pip, Pin > 0. The absolute value 
can then be replaced by \Pi\ = pip + Pin- Equivalence is 
guaranteed by minimizing the (positive) sum, such that at 
least one variable of each pair {Pip, Pin) will be zero 1(121 . 
The linearization of the cost function inherits a large error 
in the value of J, but the values of the coefficients P are 
not affected that much: the least-distance problem is not so 
much different in the linear form as it would have been in 
the quadratical form. Furthermore, see that a difference only 
appears if a constraint is active, the unconstrained optimum 
is the same. Suboptimality is evaluated in the results section 
based on numerical simulations. The linear constraints on a 
are as well transformed directly with (|T6l . Therefore, after all 
the transformations, the problem is available in standard form 
for linear programming, and a simplex solver [flSl can be run. 
The optimal solution (3* has to be re transformed to find the 
optimum in the original coordinates 

a* ^ OLo + A-^ (3* , (20) 

which is used in (fTTI) and the motor model to obtain the 
optimal trajectories for the currents and voltages. 

B. Implementation Aspects 

All computations of the previous section are done using a 
computer algebra tool (Maplesoft Maple). The matrices for the 
LP solver are generated using the C code generation toolbox. 
The real-time software thus consists of a simplex tableau 
assignment, consisting of the initialization of one matrix and 
two vectors of floating-point variables, which is automatically 
generated code, a simplex LP solver from [,18|, and some 
post-processing. As the assignment is based on symbolic 
calculations, the motor parameters can be changed online. 
However, the majority of parameters is defined as constant, 
the compiler can then optimize code size as well as runtime 
of the initialization stage. 

Alternatives to this procedure were evaluated. The direct 
use of a quadratic program is too slow, for instance, with 
6 parameters, 16 (voltage) constraints and 20 iterations it 
requires 730 /xs - if additionally current constraints were 
included, it would take even longer. If the problem is first 
transformed to least-distance, the number of iterations and the 
runtime reduce to about 50%. 

Also online calculation of the trajectory generation algo- 
rithm instead of computer algebra generated code was ana- 
lyzed. The matrix inversion with 6 parameters takes 2.5 /its and 
the cholesky decomposition 40 /is. Additionally this requires 
the use of a scientific mathematics library. 

The choice to use a computer algebra tool and automatic 
code generation turns out as the easiest way to implement and 



as well as computationally most efficient. The automatically 
generated code is copy-and-pasted to the real-time software, 
only a small LP solver like ilSlI has to be added. The chosen 
LP implementation, however, seems to be a somewhat slow 
implementation, runtime improvements are possible. 

C. Predictive Control 

The control structure is shown in Fig. |5] A cascaded 
control structure is chosen as speed is assumed constant for 
trajectory generation. As the mechanical plant is generally 
only roughly known, this overlying speed controller is advan- 
tageous. Model-based control can be used for the electrical 
subsystem of the motor as the parameters are known, but for 
the mechanical part, any robust feedback controller can be 
chosen. 

The trajectory generation scheme is embedded in a predic- 
tive controller and repeated at every sampling step to generate 
the optimal control input Udq[k]. 

Additionally, for steady-state accuracy, a disturbance ob- 
server is necessary, for instance to compensate modeling errors 
of the induced voltage terms. Even though the induced voltage 
is modeled, variations in K caused for instance by heating 
or identification errors cause an offset of the current iq. The 
approach is to assume the disturbance as constant and calculate 
it directly based on the MPC model and past control inputs 
and current measurements [[20]. The respective delays have to 
be accounted in the design. 

D. Timing 

First, measurement and control timing is analyzed. Two 
timing sequences are shown in Fig.|6](a) and (b). The interrupt- 
based control system triggers an interrupt every 125 /xs (signal 
3). At this instant, the applied voltage command (signal 2) 
is modulated in the power electronics. At the same instant, 
the A/D conversion of the current measurements (signal 1) is 
performed to avoid the impact of current ripples. It is seen that 
to calculate the command only the current iq[k — 1] is 

available. Furthermore, once the command Uq[k] is calculated, 
it needs to wait until the next interrupt. This delay of one step 
is denoted computational delay and is accounted with a delay 
compensation technique. A prediction of one sampling step 
with the model, the current iq[k — 1] as well as the previously 
commanded voltage Uq [k — 1] is applied to generate the delay- 
compensated current iq[k\k — 1] as initial condition for the 
control law to compute Uq[fc], see 1211 . Furthermore, see that 
the response to the commanded voltage Uq[k] is iq[k + 1] 
which is available one interrupt later, this is the plant delay 
and it is naturally included in the predictive controller by 
recalculating the trajectory at every sampling step. From the 
predicted trajectory, the applied voltage is Udq[k] ~ Udq{0). 

From the timing sequences, interesting insight into the com- 
putational demands of the algorithm is gained. The first part 
of the controller signal in Fig. |6]shows the calculation time for 
the simplex tableau initialization, it takes about 10/is. Included 
in these calculations, which is automatically generated code 
resulting from symbolical calculations, is a calculation of the 
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Fig. 4. TransfoiTnation and lineaiization of the optimization problem. Lines are equimagnitude contours of the cost function. Left: oiiginal problem, middle: 
transformed least-distance problem, right: transformed and linearized problem. 
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Fig. 5. Control structure of the predictive torque controller cascaded by PI speed control. 



unconstrained optimum and the linearization of the problem. 
The second and biggest part of the controller signal is the 
runtime of the linear program [18 j. At the beginning of Fig. |6] 
(a), where voltage and current are both zero, it is only about 
20yus, but to calculate the voltage step at 2000 rpm shown in 
Fig. |6] (b), more iterations are involved as many constraints 
are active, and the computation time rises to almost 60/is. 
The total time of the interrupt handling, latency, the simplex 
initialization, the LP solver and the post-processing sum up to 
almost 100/is in the worst case, therefore up to 80% of the 
available time is used. 

IV. Results 

A surface-mounted PMSM with parameters shown in table 
II is used. It is coupled to a load drive such that an arbitrary 
speed can be imposed. The algorithm is implemented on a PC- 
104 based real-time system with a 1.4 GHz CPU described 
in II22I . The voltage limitation was set to 250 V because of 
hardware limitations, accordingly, rated speed is reduced from 
3000 to 2200 rpm. 

Experimental results of the proposed scheme are shown in 
Fig. I2I Subfigure (a) shows the response to two subsequent 



speed reference steps, the load drive is deactivated. The torque 
is increased rapidly by a high but feasible voltage peak - the 
PI speed controller can be made very fast such that also a 
step appears on the reference torque for the MPC controller 
The direct current depends on the speed and thereby reduces 
iron losses which are considerable at high speeds. Losses are 
decreased by about 4%, and the efficiency is improved by 
about 0.5% at 2000 rpm. Better results are obtained on motors 
with higher inductances lITSl . 

The next three subfigures (b), (c) and (d) show fast torque 
transients at zero, medium and high speed, respectively. The 
PMSM is in torque control mode while the load drive keeps 
speed constant at 0, 2000 or 2400 rpm, respectively. The 
current components are well decoupled, a fast current change 
on the quadrature axis does not affect the direct axis in (b) 
and (c). With two separate PI controllers, a short current 
excursion would be seen on the direct axis during the torque 
transient. Again, the current on the direct axis id is dependent 
on the speed. Furthermore, the torque change is fast and at 
the same time smooth, the voltage becomes smoothly small 
for smaller control errors - a nice characteristic of quadratic 
cost functionals, compared to linear cost functions which result 



k-3 k-2 k-1 k k+1 k+2 k+3 k+4 

samples (At=0.125 ms) 
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samples (At=0.125 ms) 



(a) Timing at rpm. (b) Timing at 2000 rpm. 

Fig. 6. Experimental timing results of the predictive control scheme. 1: quadrature current iq without delay compensation, 2: quadrature voltage Uq, 3: 
interrupt handling, 4: control law computation (first part: initialization, second part: LP solver, third part: post-processing). 



in deadbeat behavior and are more sensitive to uncertainties 

lED. 

On subfigures (b) and (c), the behavior with active vohage 
constraint is the same as when using standard saturation or 
anti-windup strategies. On subfigure (d), however, a different 
behavior is seen, the direct current id is reduced to perform 
field-weakening. This implies that the stator induced voltage is 
reduced on the quadrature axis, see eq. (|3). This mechanism is 
known from minimum-time current control ll23l . With field- 
weakening, the gap between induced and maximum voltage 
increases, the derivative of the quadrature current j^iq is 
higher and the torque-generation dynamics are increased, at 
the cost of higher copper losses on the direct axis. Without 
additional field-weakening, the reference torque would not be 
reached after the optimization horizon of 2 ms, thereby the 
end-weight of the control error in J oversizes the loss term. 
Therefore, in this predictive control implementation, field- 
weakening not only improves efficiency, but also improves 
dynamics by exploiting cross-coupling between the orthogonal 
current components to optimally bypass the voltage saturation. 

The situation described above is the one where the highest 
number of constraints is active, and thus, also where the 
suboptimality of the proposed trajectory generation method 
is the highest. Based on the simulation results in Fig. |8] the 
suboptimality can be analyzed by comparing the results to 
a QP solver, which is not practicable in the experimental 
setup. In Fig. [8] (a) no real difference can be seen as only 
few constraints are active. In Fig. |8] (b) more constraints are 
active, as result, it is seen that the dynamical field-weakening 
is much more intense, but foremost, of longer duration when 
using QP. The simulations in Fig. [8] have a slightly different 
quantitative behavior as the experiments, caused by parameter 
uncertainties. 

It is also possible to operate the PMSM beyond rated speed 
with steady field-weakening to bypass the voltage saturation 
on the quadrature axis, as shown in I?]. It is remarked that 
the current on the direct axis id has no reference, its value 



is obtained from the optimization of the cost functional. 
Therefore the method works well and is numerically stable; 
the optimal value follows inherently. 

V. Conclusion 

Two questions have been central to this work, first, is it 
possible to implement linear model predictive control (MPC) 
with online optimization on an electrical drive, and second, 
will this advanced controller lead to any merits. 

A model predictive control scheme for a PMSM was 
introduced. Based on suboptimal real-time optimization, the 
currents and voltages are computed according to a cost func- 
tional. The prediction horizon is 2 ms at a sampling rate of 
8 kHz, and voltage and current constraints are respected. The 
work can be seen as first implementation of classical model 
predictive control on an electrical drive, where the aspects of 
constraint inclusion and high prediction horizon are satisfied 
while online optimization is used. Implementing online MPC 
is therefore possible even on fast-sampling systems such as 
electrical drives. 

Additionally to the advantages of general predictive con- 
trollers, which are precise accounting for timing of measure- 
ment and control as well as respecting current and voltage 
constraints, the advantages of the long-range constrained pre- 
dictive MIMO control scheme can be concluded as follows: 
improved accounting of cross-coupling, fast and smooth dy- 
namical behavior, improved power efficiency by field weak- 
ening, and improved dynamics close to voltage saturation by 
field- weakening. 

As negative point, the high demands to computational power 
must be named. 
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Fig. 7. Experimental results of the predictive control scheme. 1: rotor speed ujm, 2: direct current without compensation, 3: quadrature ciurent iq without 
compensation, 4: quadrature voltage Uq. 
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Fig. 8. Simulation results of the predictive control scheme, black: QP, light grey: LP. 1: rotor speed ojjv/i 2: direct cuiTent without compensation, 3: 
quadrature current iq without compensation, 4: quadrature voltage Uq. 
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Appendix: Motor Parameters 



TABLE II 

Nominal Parameters of the Synchronous Motor 



Manufacturer & Model 


Merkes MT5 1050 


Rated Power P]v 


2760 W 


Rated Torque tjv/at 


10.5 Nm 


Rated Cun'ent (peak) 


8 A 


Rated Speed ujmn 


3000 rpm 


Pole Pairs rtp 


3 


Rated Voltage Um (peak) 


560 V resp. 330 V 


Stator Inductance Ld, Lq 


4.8, 7.2 mH 


Stator Resistance R 


0.92 n 


Motor Constant K (peak) 


0.334 Vs 


Iron Loss Constant (Hysteresis) kp^ 


1-27 A 



